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Abstract: We investigate the effect of axions on the polarization of elec- 
tromagnetic waves as they propagate through astronomical distances. We 
analyze the change in the dispersion of the electromagnetic wave due to its 
mixing with axions. We find that this leads to a shift in polarization and 
turns out to be the dominant effect for a wide range of frequencies. We 
analyze whether this effect or the decay of photons into axions can explain 
the large scale anisotropies which have been observed in the polarizations 
of quasars and radio galaxies. We also comment on the possibility that the 
axion-photon mixing can explain the dimming of distant supernovae. 

1 Introduction 

The mixing of axions with photons and its observable consequences have 
been analyzed by many authors [1-12]. In the present paper we investigate 
the changes in the polarization of electromagnetic waves that arise due to 
its mixing with axions. We are particularly interested in determining if this 
mixing can explain the polarization anisotropies that have been claimed in 



Ref. 13, M 



In Ref. [D5, [nj the authors claimed that the observed polarizations from 
distant radio galaxies and quasars are not isotropically distributed on the 
dome of the sky. The observable of interest in that study was the angle 
(3 = x ~ where ip is the orientation angle of the axis of the radio galaxy 
and x the observed polarization angle after the effect of Faraday rotation is 
taken out of the data. The authors claimed a dipole anisotropy such that 
the angle j3 is given by 

p = X-r (1) 
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where f is a unit vector in the direction of the source. The A represent the 
three parameters of this fit. The magnitude of this vector |A| is found to be 
approximately 0.5 and its direction 



A = [(0 h,9m)± (1 h, m), -1° ± 15°] 



(2) 



The effect is independent of redshift and was first claimed in Ref. |13| and 
later verified by a more reliable statistical procedure |T7| and by compiling 
a larger data set fI3| . This anisotropy may be a signal of some local effect 
arising due to the milky way or the local supercluster. However so far it is 
not known what physical phenomenon could lead to the observed rotation 
in polarizations. Within the standard model of elementary particles it is 
difficult to conceive of a physical mechanism which can lead to this effect. 
The axion field which arises in many extensions of the standard model of 
particle physics may provide one possible explanation. However so far it is 
not known whether such a field can consistently explain this effect. This is 
one of the questions that we study in the present paper. 

Another interesting polarization effect in the electromagnetic waves from 
distant quasars has been claimed in Ref. [Tj| IB|. It was found that opti- 
cal polarization are aligned on very large scales. A very striking alignment 
was found in the region, called Al in jTjJ, delimited in Right Ascension by 
ll h 15 m < RA < 14 h 29 m and in redshift by 1.0 < z < 2.3. The polarizations 
from quasars in any particular spatial region have a tendency to align with 
one another. The effect was only seen in patches without any evidence of 
large scale anisotropy. We point out that the center of the Al region (see 
Fig. 1 of Ref. |16[) is exactly opposite to the axis, Eq. 0, of the anisotropy 
found in Ref. |14| . This might indicate a common origin of these two effects. 



Finally we also examine the recent claim that dimming of distant 
supernovae |19|, |20[ can be explained in terms of axion photon mixing. 



2 Axion-Photon Mixing 

The interaction lagrangian of the axions with electromagnetic field can be 



written as 21 



c 



int 



aN 



(3) 



where is the axion field, F„ v is the electromagnetic field tensor, f a is the 



scale of PQ [22] symmetry breaking and N is the number of light quark 
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flavors. The current limits on f a are given by 




(4) 



(5) 



where f n is the pion decay constant, m n the pion mass, and m u and the 
masses of up and down quarks respectively. It is also interesting to consider 
other pseudoscalar particles which arise in certain extensions of standard 
model, whose mass is not related to f a by Eq. |]. In most of our discussion 
below we will take the mass of the pseudoscalar as a free parameter and not 
given by Eq. [5|. 

Axions can mix with photons from distant galaxies and lead to rotation 
of polarization. The basic picture is that the photons emitted by the galaxy 
can decay into axions as they propagate through the background magnetic 
field. Since only the photons polarized parallel to the transverse component 
of the background magnetic field (B t ) decay, this effect can lead to a change 
in polarization of the electromagnetic wave. Another effect that can also 
contribute is the mixing of photon with an off-shell axion. This process 
changes the dispersion relation for the photon polarized parallel to B t and 
hence can lead to changes in polarization. Alternatively the distant galaxies 
may be emitting axions along with photons. These axions can decay in the 
presence of magnetic fields into photons which are polarized parallel to B t . 
We examine all of these possibilities. 

The basic equations can be written as || || , 



where u p is the plasma frequency, A\\ is the component of the vector potential 




(6) 




(7) 



parallel B t , B, 



•t — 



B t \ and the coupling g 



aN 
'in f a 



. By making the ansatz, 
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and assuming that the frequency u is much larger than the mass eigenvalues, 
the wave equations can be written as, 



id z A\ + ^ exp(-i^(z) /2u)(f) / = 



id z (f>' + =^exp(i£(z)/2a;)Aj| = 



where 



dz\u? v {z') - ml) 



(10) 

(11) 
(12) 



In arriving at equations [K| and |TT] we have ignored second derivatives of A'n 
and 4>' since these are slowly varying fields. 

We first replace the ujp(z) with its mean value Up. In this case, as shown 
in Ref. |J, the probability for producing the pseudoscalar particles |0(L)| 2 
at distance L, assuming that 0(0) = is equal to 



|0(L)| 2 = |0'(L)| 2 ^(^/) 2 sin 2 



L 
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where 



2oo 



out 



mi 



(13) 



(14) 



This gives a very small value assuming the current limit for the value of 
g and the galactic values for B t and U p , independent of the choice of m a . 
The special case where a; 2 is accidently equal to m 2 is not considered in this 
paper. The effect is, however, much larger if the variations in the plasma 
frequency and/or background magnetic field are taken into account 0. The 
authors in Ref. assume a Kolmogorov power spectrum, P 3 n = Cffk~ u / 3 
for the electron density fluctuations, where Cn ~ 3 x 10 -4 m _20//3 for the 
interstellar medium. The authors then show that for the electromagnetic 
wave propagating through interstellar medium the probability to produce 
is given by, 



axions, P 7 



11/3 



2.7 x 10 



x 




-11/3 



(15) 
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In obtaining this result the magnetic field has been assumed to be constant. 
We have also assumed that the vector potential is approximately constant 
and kept only the leading order term in the expansion in powers of the fluc- 
tuations. Hence this result is valid only if the conversion probability is small. 
The conversion probability given by this equation is much larger than that 
implied by Eq. [H^. For propagation over galactic distances this probability is 
still very small compared to unity unless the frequency of the electromagnetic 
wave is much larger than the optical frequencies. However for supercluster 
magnetic fields the probability may be equal to unity even at optical fre- 
quencies. For example, in the Virgo supercluster of galaxies the magnetic 
field strength is found to be about 1 fiG over a very large length scale of 
10 Mpc with plasma density of the order of 10 -6 cm" 3 |23|| . In order to 
compute the correlation coefficient Cn for the supercluster plasma density 
we assume, by dimensional analysis, that it scales as nf/ 9 , since n e is the 
only dimensionful parameter known. This gives Cfj ~ 3.4 x 10~ 14 m _20//3 . 
This estimate of is not reliable and must be improved in future by di- 
rect observations. Using this rough estimate the probability is found to be 
approximately 8 x 10 5 for g = 10~ 10 GeV -1 . This is large enough that it 
can also significantly affect the CMBR in the direction of the Virgo super- 
cluster. By demanding that < 10~ 5 for the microwave photons v = 50 
GHz we find that g < 0.9 x 10~ n GeV -1 . This constraint is comparable to 
the constraint on coupling g obtained from red giants and SN1987A [0, ||. || 
but depends on our assumed extrapolation of Cn and hence is not entirely 
reliable. The constraint imposed by CMBR also implies that 

< 6 5 X 10 (l^Gifa) (16) 
The conversion probability is of course always less than one. The large con- 
version probability is obtained by using Eq. [15] in the regime where it is not 



valid. As mentioned above, Eq. [15] is applicable only when the conversion 
probability is small. In the present case the result, Eq. [16], has to interpreted 
in the sense that for a considerable range of parameters P 7 _»^ is of order unity 
for v > 10 6 GHz. 

We next investigate how the polarization of photons will change due to 
change in the dispersion of photon because of their mixing with axions. We 
integrate Eq. [ll] by parts, taking B t as constant, to obtain 



- 0'(O) = -i^-l 



* /l AUz) - AUO) - f dz'd z A\ exp(iz'/l) 
ii ii j Q ii 



(17) 
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It is clear that the third term inside the brackets is higher order in g and can 
be dropped. Substituting this expression for 0' into Eq. [10] we find that 



d z A' n 



-i^Y<j>'{0)exp(-iz/l) +i 



gB t 



AUz) - ex-p(-iz/l)AU0) 



(18) 



The solution to this equation can be written as, 



x 



Aj|(0)exp(zr(L)) - 

2^(4 + ^ 
\2l 2 



I 



1 + W(L)/L 



gB t 



0'(O) + ^L4(,(O) 



+ i sin 



exp(zr(L)) 



where, 



r(L) 



(gB t 



V 2 



IL 



(19) 



(20) 



For large L, of the order of galactic or super-galactic scales, the dominant 
contribution comes from the first term on the right hand side of Eq. This 
term leads to an additional phase T for the parallel component of the vector 
potential in comparison to the perpendicular component. The total photon 
flux fluctuates due to its conversion into axions. 

In the optical frequencies, assuming the galactic values of the magnetic 
field B t ~ lfiG and the plasma density n e « 0.03 cm -3 , we find that T(L) 
is of order of 10 -5 for L w 1 kpc and g w 10~ 10 GeV -1 which is comparable 
to the probability of conversion of photon into an axion as given by Eq. [1"5| . 
However for smaller frequencies this gives a much larger result in comparison 
to Eq. [T5|. If we assume the Virgo supercluster values of the magnetic 
field B t ~ lfiG and plasma density n e w 10 -6 cm" 
m a « uJ p , 

0.04 Mpc 



10 6 GHz 

and the phase T can be written as 



| 23| we find that for 
(21) 



r(L) 



2.4 x 10' 



g 



10- 10 GeV- 




x 



L 


V 


V 10 Mpc/ 


L10 6 GHz J 



(22) 



This phase can produce observable changes in the polarization of the elec- 
tromagnetic wave as we discuss later. 
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We next compute the contribution to the phase that arises due to fluctu- 
ations in the plasma density. We assume, for simplicity, that the background 
magnetic field is constant. Then solving for <fr' in Eq. [H] and substituting 
into Eq. [10] we find, 

Mi = - (t0 2 -<> (-4) /» Vexp (4) A > (23) 

Here we have set <p'{0) = 0. If <p'(0) ^ then we will find another term 
on the right hand side which can be evaluated easily to leading order in the 
fluctuations. We discuss this case later. We next write A'u = e in ~ 1 + zO, 
where Q is in general complex. We set -Am(O) = 1 and hence f2(0) = 0. 
We assume that Q(L) is small and keep only the leading order terms. In 
order to compute the higher order terms we will require the higher order 
density correlations functions for the turbulent interstellar (or intergalactic) 
medium. These correlation functions are, however, unknown and hence it 
is not possible to go beyond the leading order term. In any case even the 
leading order calculation is very useful since it can reliably provide us with the 
parameter ranges where the axion photon mixing effect may be significant. 
We find that Q(L) can be written as 

aiL) = ' (*r )7„ Wp (-4) (4) < 24 > 

We can now expand £/2u; in terms of the fluctuations in the plasma density 
to get 

f = T + f- f d^ne (25) 
Zuj i Zuj n e Jo 

where we have used ioi oc n e , = Td 2 v + 5oj 2 and n e = n e + 5n e . Substituting 
this into Eq. |24] and expanding the right hand side to second order in the 
fluctuations 5n e we can write the expectation value < Q >= Qq + ^ with 

where, 

I\ = lj dze %z/l 5n e J dz x bn e 
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I 2 = ^jf dze" iz ' l 5n e £ dz ie izi/l 5n P 
h = ^ dze- iz/l 5n e J* d Zl 5n}j (27) 

and Qq is the contribution obtained by keeping only the mean electron density 
and can be extracted from Eq. |19[ 

These integrals can be evaluated by expressing them in terms of the 
Fourier transform of 5n e = J c/ 3 kfi e (k)e* k ' x and using the Kolmogorov power 
spectrum (n e (k)n*(k')) = P 3 jv(£;)<5 3 (k — k'), where k = |k| and 

P Mk) = (fc2 + ^ 2)11/6 exp (-kHl/2) . (28) 

Here is the correlation coefficient, Iq is the inner scale and L the outer 
scale f24|. In our calculation we set l = 0. We have explicitly kept the 
dependence on the outer scale Lo since we find that some of the integrals are 
infrared divergent which are regulated by this scale. We point out that since 
5n e is real n*(k) = n e (— k). We find that 5fl is given by 



oo 




L - \< > > (29) 



In obtaining this result we have kept only the dominant terms which scale 
like L in the limit of large L » l,L » L Q . The imaginary part Qi of 
< Q >= Qfi + iQj, gives the decay rate of the photons with polarization 
vectors parallel to the transverse component of background magnetic field. 
It is related to the production rate of axions by \<f)'(L)\ 2 = 2f2/. The real 
part, Qr, gives an extra phase factor to the parallel component. The integral 
in the real part of 5Q can be done numerically. We find, 

It is clear from Eq. [30] that the Sftji/Sftj ~ (L //) 2//3 , where 5Qr and 
5Qj are defined by 5Q = 5Qr + iSQi. Hence we find that SQr is much larger 
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than Sfli if the frequency is such that I << L, I << L and L << L. 
We therefore find that for a wide range of frequencies which satisfy these 
inequalities the phase factor Qr » Qi irrespective of whether the plasma 
density is uniform or has a fluctuating part. In both cases we find that the 
contribution to Qr qc uj. For the Virgo supercluster, assuming that Cfj scales 
as n e and L scales as L we find that 5Qr ~ 9 x 10 5 for optical frequencies 
v = 10 7 GHz and g = 10~ 10 GeV -1 . Hence at optical frequencies the real 
part of 5Q is comparable to the imaginary part. However at lower frequencies 
the real part is much larger. 

We next consider the case when the source is emitting axions along with 
photons. In this case <f)'{0) 7^ 0. The expression for (L) can now be written 
as 

A\(L) = -i^'(0) £ dzex P Li£\ 

- £* fe *>{ i l) (31) 

In order to compute the polarization properties of the electromagnetic wave 
after propagation through a large distance L we need to compute the corre- 
lations < A'u(L)AUL) > and < A±(L)A'^(L) >. We ignore all terms which 
involve third and higher powers of the coupling g. Furthermore we include 
only those terms which scale like L for very large values of L » I. With these 
requirements it is clear that only three terms contribute to < A'u(L)A'JL) >. 
These are |Aj,(0)| 2 , 




^ (33) 

The second term (Eq. ^) given above is equal to the result given in Eq. 
upto an overall factor of |0'(O)| 2 . The third term, given in Eq. |33], is equal 
to the RHS of Eq. ^4] (plus its complex conjugate) up to an overall factor 
of zyljf (0)ylj|(0). The fluctuating part of this expression is evaluated in Eq. 
j30[ . Similarly the only nontrivial term that contributes to < A^(L)A',T(L) > 
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involves the expectation value of the third term on the RHS of Eq. [JT] which 
is already evaluated in Eq. |3(J We point out that the expectation value of 
the second term on the RHS of Eq. does not scale like L in the limit of 
large L. Hence the terms which involve interference of this term with a term 
independent of fluctuations is neglected. We, therefore, find that the only 
change introduced due to 0'(O) 7^ is an extra term in < A!»(L)AUL) > 
which arises due to decay of axions. This is the term given in Eq. |32| and 
upto an overall factor of |0'(O)| 2 , its expression is given in Eq. [15]. 

We next compute the changes in polarization induced by the phase Q R 
due to axion photon mixing. We neglect the imaginary part of Q in the 
following discussion. As shown above it is in general much smaller than the 
real part. Furthermore its effect can be easily included by multiplying the 
component of vector potential (or the electric field) parallel to B t by e~ n ' . 
If the incident electromagnetic wave is unpolarized then the phase fl R does 
not produce any additional effect. However if the incident beam is already 
polarized then Q R can change also its polarization. Let's consider the simple 
case of a monochromatic wave, 

E x = ai exp(i«i — iuot) 

E y = a 2 exp(ia 2 — iujt) (34) 

where E x and E y are the two transverse components of the electromagnetic 
wave. By choosing the coordinate axis such that the y-axis is aligned along 
the direction of the transverse component of the background magnetic field, 
we find that after propagation through the medium the y component of the 
electric field E y acquires an additional phase Q R . Hence the final values of 
the Stokes parameters are given by 

50 = a\ + aj (35) 

51 = a\ — a\ (36) 

5 2 = 2a 1 a 2 cos(A + Q R ) (37) 

5 3 = 2a x a 2 sin(A + Q R ) (38) 

where A = a 2 — The linear polarization angle ip is then given by 

S 2 2a l a 2 cos(A + n R ) 

tan 2</> = — = - 2 — - 2 (39) 

01 a x — a 2 

The degree of polarization P = (Sf + S% + iSf) 1 ' 2 /^ remains equal to unity. 
Similarly the degree of circular polarization defined as Pc = \S$\/Sq is given 
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by 

P c = \2a 1 a 2 sm(A + n R )\/(a 2 1 + a 2 2 ) (40) 
and the degree of linear polarization Pl = (Sf + S 2 ) 1 ' 2 / Sq is given by 

P L = {(a 2 - a\) 2 + Aa\a 2 2 cos 2 (A + Q R )]^ 2 /3 • (41) 

We see that the linear polarization angle ip, the degree of circular polarization 
Pc and the degree of linear polarization Pl all depend in a precise manner 
on the frequency uj of the incident wave. If we assume that Qr is small then 
these depend linearly on uj. In general these are oscillatory functions of uj. 
On the other hand the degree of polarization P is independent of frequency. 
Hence we can easily check this effect by analyzing spectral data from distant 
AGNs at optical and higher frequencies. 

If the incident beam is quasi-monochromatic the above analysis goes 
through with minimal modification. The wave now need not be completely 
polarized. We may write its coherency matrix as 



Jn J12 

J 21 J22 



(42) 



where Jij =< E*(t)Ej(t) > and J 2 i = J*2 P5| - After propagating through 
the magnetic field which is assumed to be pointing in the y direction we find 
that the y-component E 2 of the electric field acquires an additional phase 
Qr. Hence the coherency matrix changes to 



12 

21 ^22 



r m «J a1 J, 



(43) 



The Stokes parameters can be written as So = Jn + ^22 > Si = Jn — J22, 
S 2 = e tQR Ji2 + c.c. and S 3 = i(e~ inR J 21 — c.c). Parametrizing J 12 = je lS we 
find that S2 = 2j cos(5 + Qr) and S3 = 2j sin(5 + Qr). We again find that 
the degree of polarization is independent of Q but the orientation of linear 
polarization as well as the degree of linear and circular polarization (Pl and 
Pc) become oscillatory functions of u due to their dependence on Q R . 



3 Astrophysical Applications 

We next examine if the astrophysical effects discussed in the introduction can 
be explained in terms of axion photon mixing. As we have seen the effect 
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is too small for the interstellar magnetic fields. However the effect can be 
large and produce observable consequences for the case of Virgo supercluster 
magnetic fields, especially at optical frequencies. We first consider the po- 
larization alignment effect claimed by Hutsemekers |15|. This effect can be 



explained by the depletion of the photons polarized parallel to B t . We point 
out that the real part of Q by itself cannot explain this effect. As explained 
in the introduction, a very striking alignment was found in the region Al in 
ref. [113], which is bounded in Right Ascension by ll h 15 m < RA < 14 h 29 m 



and in redshift byl.0<z<2.3. In Ref. |L6|] it is also shown that the quasar 
polarization in the region Al show alignment with the supergalactic plane. 
By invoking coherent cluster magnetic field it may be possible to explain the 
existence of large scale correlation. The region Al turns out to be in the 
direction of the Virgo supercluster and this appears to be a promising ex- 
planation. As pointed out in Ref. |T?J the correlation with the supergalactic 
plane may be evidence in favor of a propagation effect such as axion photon 



mixing. An alternate possibility, also mentioned in Ref. |]16||, is extinction 



due to dust. However this explanation suffers from two drawbacks. The first 
problem is that the effect is redshift dependent, i.e. the regions which show 
alignment are bounded in redshift as well as in angular coordinates. Sec- 
ondly one finds that in general the distribution of the degree of polarization 
for the radio quiet (RQ) along with optically selected non-broad absorption 
line quasars (O) differs from the distribution of the broad absorption line 
(BAL) quasars. Precisely the same difference between these different types 



of quasars is also seen in the Al region [16]. If the supercluster magnetic field 



causes the alignment in polarization due to decay of An then this difference 
would not be preserved in the Al region. 

It, therefore, seems that a propagation effect is not likely to provide an 
explanation for this effect. This would, however, imply that the quasar polar- 
izations are intrinsically aligned with one another over cosmologically large 
distances. Correlations over cosmological distances violate the basic assump- 
tion of isotropy and homogeniety of the Universe. Here we propose an alter- 
nate mechanism which explains these observations in terms of a propagation 
effect related to axion photon mixing. We first notice that there are very few 
RQ+O quasars in the Al region which satisfy the cut p > 0.6 % imposed 
on the data [|15|. Hence these quasars do not contribute significantly to the 
alignment effect. We next assume that the quasars are emitting axions along 
with the photons. During propagation the axions decay into photons in the 
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presence of background magnetic field. The decay probability of axions is 
also given by Eq. [15]. The emitted photon is polarized in the direction par- 
allel to the transverse magnetic field. The axion photon coupling is assumed 
to be such that in the optical frequencies the decay probability of axions, 
or that of photons, from distant quasars is small, less than about 1 %. The 
total photon flux generated due to decay of axions is equal to the product of 
the initial axion flux $ and its decay probability P^j. We further assume 
that for the BAL quasars the initial axion flux $ is large enough at opti- 
cal frequencies such that it leads to significant change in polarization during 
propagation. In order to explain the Hutsemekers effect we require that the 
axion flux is of the same order of magnitude as the photon flux from BAL 
quasars, if the probability P^ 7 ~ 1 %. In the current paper we will not 
address the question of emission rate of axions from quasars since the basic 
physics of the interior of the quasars is not well understood. For the (RQ+O) 
quasars we assume that the emitted axion flux is small enough such that it 
does not affect the optical polarizations. The distribution of the degree of 
polarization for the (RQ+O) quasars will, therefore, be unaffected due to ax- 
ion photon mixing. However the degree of polarization for the case of BAL 
quasars will be determined primarily by the propagation effects. Hence this 
mechanism explains the observed difference between the BAL and (RQ+O) 
quasars. 



The redshift dependence of the Hutsemekers effect [15, 16 1 is difficult to 



explain even with the axion decay mechanism discussed above. However 
we notice that the dominant contribution to the effect comes from the Al 
region which can be explained in terms of the Virgo supercluster magnetic 
fields. The strong alignment of the polarizations with the supercluster plane, 
as shown in Fig. 3 of Ref. |fL6"| , provides further support for this explana- 
tion. Some redshift dependence may be obtained due to the evolution of the 
quasars with redshift. It is possible that the axion emission rate of quasars 
is redshift dependent and hence only quasars within a certain redshift in- 
terval give dominant contribution, with the remaining quasars dominantly 
contributing noise. It is also possible that the dominant contribution to 
the effect comes from a few clusters within our astrophysical neighbourhood 
which can then also contribute to the observed redshift dependence. This 
can be studied by collecting a larger data sample and putting appropriate 
cuts on the redshift. 

It is interesting to note that the observed difference between the BAL and 
RQ+O quasars in the direction of the Virgo supercluster allows us to put a 
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very stringent constraint on the pseudoscalar photon coupling g in the limit 
of m a << ZJ p provided we assume our extrapolated value of the correlation 
coefficient Cn- We may isolate the influence of the Virgo supercluster by 
taking only the region in which a large Faraday Rotation Measure is observed. 
This corresponds to a patch of 20° radius around the center of the Virgo 
supercluster (RA=12 h 28 m , DEC=12°40') |23|. We find that in this region 
there are 16 BAL quasars and 9 RQ+O quasars among the objects given 
in table 2 of Ref. [03]. The difference between these two types of quasars 



persists in this region also as can be seen in Fig. |T]. If the photon to axion 
conversion probability is large, i.e. P 7 ^ > 0.1, then this difference would be 
completely washed out due to propagation through the Virgo supercluster. 
This is because the degree of polarization p of RQ+O quasars peaks at very 
small values (p ~ 0.2 %) and we find no quasars in this category with p > 2 
%. This observation implies that the photon decay probability P- 



< 1, 



where the upper limit, unity, is chosen in order to obtain a conservative upper 



bound on g. By using Eq. 15 and the parameters for the Virgo supercluster, 



given in the discussion following Eq. 15, we find that 



g < 10~ 13 GeV 



3.4 x 10 



-ii 



m 



-20/3 



r<2 



1 fiG 



10- 6 



cm 



10 Mpc 



L 



(44) 

The reliability of this constraint depends on our assumed extrapolation of 
the correlation coefficient Cjy for the supercluster. In future it is clearly of 
interest to obtain this parameter by direct observations which can lead to a 
more reliable constraint. 

In the radio frequencies (~ 1 GHz) the effect is clearly very small. The 
dominant contribution is given by 5Qr assuming Virgo supercluster param- 
eters, given in the discussion following Eq. We find that for radio fre- 
quencies v = 1 GHz, 5Vt R « 0.1 for g = 10~ 10 GeV" 1 with m a « uJ p . The 
appropriate value of g is, however, likely to be much smaller than g = 10~ 10 
GeV -1 . We next investigate whether the proposed emission of axions due 
to quasars can explain this effect. Given a large enough flux of axions it is 
clearly possible to obtain an observable effect on the polarization at radio 
frequencies. The probability for decay of axions into photons P<£_> 7 in the 
presence of background magnetic fields is also given by the result given in 
Eq. [15|. Using the CMBR constraint on the coupling g we find that the 
P(j>-n — 6-5 x 10 -7 at radio frequencies v ~ 1 GHz. Hence in order that 
axions decay can lead to observable consequences the axion flux has to be 
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0.4 0.8 1.2 1.6 2 

Polarization Degree 

Figure 1: The distribution of polarization degree p (in %) for the BAL quasars 
(solid line) and the RQ+O quasars (dashed line) for a patch of 20° radius 
centered around the Virgo supercluster. 
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roughly 10 times larger than the photon flux F y at radio frequencies. The 
photon flux at radio frequencies is given by Cjv q where C is a constant and 
the spectral index q varies roughly between 0.5 to 2 and is in general larger for 
larger frequencies. Assuming that the axion flux from quasars is indeed 10 6 
times larger than the photon flux at radio frequencies and that its spectral 
dependence is given by A/v p , where A and p are constants, we find p > q+5/3 
in order that photon spectrum is not in conflict with observations. Assuming 
that q ~ 0.8 at radio frequencies we find that p ~ 2.5. The axion luminosity 
of the quasar is given approximately by A7rL 2 A/[(p — l)z/^ _1 ], where v c is the 
lower cutoff on the frequency and L is the total distance of propagation. We, 
therefore, find that the axion luminosity is larger than roughly 10 6 times the 
radio photon luminosity. For most quasars this is one to two orders of mag- 
nitude larger than their total bolometric luminosity. This large luminosity 
required suggests that this mechanism is disfavored and we do not pursue it 
further. 

Another mechanism that might be relevant for explaining the anisotropy 
in radio waves involves the presence of background axion field. As shown in 
Ref. Q the rotation in the linear polarization angle due to a background 
magnetic field is equal to gA<f)/2 where A0 is the total change in the back- 
ground axion field along the trajectory of the electromagnetic wave. If we 
assume that the axion field in our neighbourhood shows a dipole distribu- 
tion, 0(f) = O (r) cos(0) where 9 is the angle with respect to the dipole 
axis and 0o(O) is of order 1/g, then the dipole anisotropy claimed in Ref. 



13, 14] can be explained. This explanation requires that the axion field at 



very large distances, which correspond to the positions of the quasars, is 
relatively smooth. The random component of the axion field should have 
a strength smaller than 1/g. However it is difficult to theoretically justify 
the existence of such a dipole distribution of the axion field and hence this 
explanation is also not very compelling. 

Finally we examine the recent claim [[18] that the observed supernova 



dimming fl9| , |20|j at large redshifts can be explained in terms of axion photon 
mixing. It has been pointed out that axion photon mixing is unlikely to 
provide an explanation because inclusion of the intergalactic plasma density 
leads to a considerable reduction of the effect Ppl. Furthermore it also gives 



rise to a spectral dependence which is incompatible with observations f26 [. 
The intergalactic parameters used in Ref. J26| are B = 10 -9 Gauss and 
n e = 10~ 7 cm -3 . However in Ref. [27| the authors argue that the intergalactic 
plasma density is likely to be much smaller than this. By using n e = 10~ 8 
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cm 



3 it is found that the axion photon mixing can provide an explanation for 

Here we examine the magnitude of the effect due 



the dimming effect [27, 28 



to the fluctuating part of the plasma density as given in Eq. [15|. Here again if 
we assume that L c << I then the I, and hence the spectral, dependence goes 
away. The result in this case can be obtained by calculating \<j)'(L c )\ 2 which is 
the probability P 7 _ + over a single coherence length and then adding all the 
contributions over the entire propagation distance. The decay probability 
over a single coherence length is given by H 



gB t u 



2n e (ul - m 2 a ) 



< I > 



where 



L c ( U 2 — m 2 ^ 

dz5n e exp i— -z 

o \ 2uj i 



We then find, using I » L c , that 
< I >= A.ix I dk. 



- - am*(k g L e /2) 

T 



dkl ~~~ v'-'x-' p.. 



3JV 



(45) 



(46) 



(47) 



Here P^n = C 2 N /{k 2 z +k^+ l/L 2 ) 11 ^ where the outer length L has been set 
equal to the coherence length. We perform the kr integral to find 
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(48) 



The remaining integral can now be done numerically and is found to be 
approximately equal to 2. It is clear that the effect is independent of v. 
Adding the contributions from the entire propagation distance we need to 
replace L^ 3 LL b J z '. Here we have assumed that the total effect is small and 
the initial axion flux is negligible. We have also ignored the fact that magnetic 
fields in different domains are aligned in different directions. A more reliable 
summation over different domains in obtained in Ref. [25|. Here we are 
only interested in an order of magnitude estimate and hence these additional 
refinements are unnecessary. In order to estimate the magnitude of the effect 



we take L = 1 Gpc, L c = 0.1 Mpc, g = l(T n GeV 



10- 



cm 



and 



again assume that Cfj scales as n 20//9 . We point out that the parameters 
L c , g and are subject to considerable uncertainty. We have chosen a 
range which is allowed by current observations and for which the effect is 
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significant and independent of v. With these choice of parameters we find 
that P 7 _>0 ~ 1.5, i.e. the conversion probability is of order unity. It is clear 
that the range of allowed parameter space is even larger if we use n e = 10 -8 
cm -3 . Hence we cannot exclude the possibility that axion photon mixing 
may be responsible for the dimming of distance supernovae. 



4 Conclusions 

In conclusion we have analyzed the mixing of axions with photons. We have 
determined how the dispersion relations of the electromagnetic waves are 
modified due to the presence of axions. We found that for a wide range 
of parameters this provides the dominant contribution to the changes in 
polarization of electromagnetic waves due to mixing with axions. We have 
also determined whether the existence of axions can explain the large scale 



anisotropies claimed in References [jn|, [14], [La]. We find that the Hutsemekers 



T5fl effect may be explained by the supercluster magnetic fields if we assume 



that quasars emit axions such that their flux is of the order of or larger than 



the photon flux at optical frequencies. The Birch effect JIB], [14]] may also 
be explain by assuming axion emission from radio galaxies and quasars but 
requires a very large flux at radio frequencies. Its explanation in terms of 
axions is therefore disfavored. Finally we have estimated the contribution of 
the fluctuations in the plasma density to the dimming of distant supernovae 
due to axion photon mixing. We find that this effect can be rather large and 
for a considerable range of currently allowed parameter space can provide an 
explanation for the supernova dimming. 
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